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This paper (lives an algorithm jrom which an approximation to the 
amplitude of a quadrature pair can be obtained simply and without re- 
quiring more bits than are needed to represent the largest allowable ampli- 
tude, A. 

If A cos 9 and, A .sin 9 are a quadrature pair, the amplitude is 

A = (A 2 cos 2 6 + A 2 sin 2 0)*. 

If AX is the larger magnitude of the pair and AY the smaller, an ap- 
proximation to A is 

A' = AX + \AY. 

Computations using this algorithm with random values of 9 showed that 
when the (rue amplitude was A, the mean obtained was 1.087 A, with stand- 
ard deviation 0.031 A. 

When A consists of a signal appreciably contaminated with noise so 
that several estimates are required before the signal can be detected reliably, 
no significant degradation («0./ dB S/X) in detectability was found using 
the algorithm described above. 

I. INTRODUCTION 

When a discrete Fourier transformation is applied to samples from 
a data segment the results are usually desired in the form of phase 
and amplitude estimates of the components into which the data is 
resolved. When FFT 1 algorithms are used to perform the transforma- 
tion the information from which these estimates are obtained is 
usually produced in the form of quadrature pairs of numbers A cos 9, 
and A sin 9, where A is the amplitude and 9 the phase of the estimate. 
Obviously 

A = (A 2 cos 2 9 + A 2 sin 2 0)*. (1) 
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The intermediate step in this computation of squaring the com- 
ponents of the pair requires in general twice as many bits as are 
required to represent A. Computing the squares requires appreciable 
time as does the final square root computation. The situation is ag- 
gravated by the fact that there are as many pair members as the 
original data samples. 

This paper describes the use of an algorithm which gives an 
approximation that in many if not all cases is adequately accurate 
and greatly reduces the computation time as well as promoting effi- 
ciency in the use of hardware by not requiring more bits than are 
needed to represent the final estimates. 

II. DISCUSSION 

Let the larger of cos 2 0, sin 2 0, be represented by X 2 and the smaller 
by Y 2 . Then 

A = AX (l + |^j 

- J *( 1 +iF-iP+--)' (2) 

where 1 ^ Y/X. 

An approximation to A is thus 

A' = AX + IAY (£) , (3) 

which suggests the computationally more attractive form 

A' = AX + \AY. (4) 

Note that AX and AY are magnitudes and (4) implies that the ampli- 
tude of a Fourier component is approximately the sum of the larger 
and half the smaller of its quadrature pair magnitudes. 

III. EVALUATION 

Curve (a) in Fig. 1 shows a plot of (4) as 6 ranges from zero to 
tt/4 where the curve is reflected to 6 = tr/2. The pattern is then re- 
peated for the remaining three quadrants. The corresponding plot for 
(3) is shown in curve (b). 

From the samples displayed on this figure the mean and standard 
deviations were calculated and found to be: 
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For (3), mean = 1.01, standard deviation = 0.0157, and 
for (4), mean = 1.087, standard deviation = 0.031. 

The algorithm described by (3) thus gives more accurate estimates 
(as expected) than that described by (4) which, however, is easier to 
implement and requires less computation. 

When the Fourier analysis is carried out in an attempt to detect 
signal components contaminated by noise it is important to determine 
how processing techniques affect the detectability of such signals. In 
order to check this a program was prepared to compute performance 
estimates when the algorithm based on (4) was used. These estimates 
were plotted on standard ROC curves 2 for a linear detector and no 
significant degradation in detectability was observed. Extensive com- 
putations using the proposed algorithm have produced spectrograms 
of noise-contaminated signals that are undistinguishable from those 
obtained conventionally. 

Use of the algorithm described in equation (4) results in a phase- 
dependent variation of the detector output. However, when the phase 
is uniformly distributed in the range (0, 2v) the amplitudes are quite 
densely clustered around the mean because of the shape of the modu- 
lation characteristic. Eighty-five percent of the amplitudes lie within 
a range of about 0.4 dB. The accuracy obtained using the fast ampli- 
tude algorithm based on (4) may thus be adequate even when the esti- 
mates are not corrupted by noise. 

IV. CONCLUSION 

An algorithm has been described for computing the approximate 
amplitudes of Fourier harmonics from their quadrature pairs. The 
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Fig. 1 — Algorithm modulation characteristics. 
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algorithm is fast and does not require the use of more bits than are 
required to represent the largest permissible result. Use of the 
algorithm does not significantly degrade the detectability of signals 
contaminated by noise, compared with use of a perfect linear detector. 
Estimates obtained using the algorithm are amplitude modulated 
depending on their phase but the total modulation range is less than 

1 dB. Representing the magnitude of the larger component of a quad- 
rature pair by AX, and the smaller by AY, the algorithm is 

A' = AX + \AY. 
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